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We analyze the supersymmetric contributions to the B — > Kit process. We show that the simulta- 
neous contributions from the penguin diagrams with chargino and gluino in the loop could lead to a 
possible solution to the B — ► Kir puzzle. Our result indicates that including the stringent constraint 
from the 6 — > sy branching ratio, the supersymmeteric models with light right-handed top-squark 
and large mixing between second and third generation of up and down squarks is the most preferred 
by the current experimental data. 



Recent experimental measurements for the CP-averaged 
branching ratios of B — > Kir decays exhibit a possible 
discrepancy from the Standard Model (SM) prediction 
[2,3]: 



R, 



Br\B A 



Br[B~ -> K-tt°] 



Br[B+ -> K°tt+ 
:i.l2±0.05) exp , 
p _ 1 (Br[B° -> K+tt- 
'" = 2 \ W[B° K"tt°} + BR[B° 

(0.79±0.11)exp- 



Br[B~ -> K°tt-] 
-Rr~(B Q -> K-TT+ 



(1) 



(2) 



As discussed in [4] , it is very difficult to have the situation 
R c > 1 and R n < 1 within the SM. While a confirmation 
with more accurate experimental data is necessary (an 
overestimate of 7r° can lead to a similar pattern of R c 
and R n , see more in detail [5]), the current experimental 
values in Eqs. (1) and (2) does not seem to be possible 
even if we consider large hadronic uncertainties [6,7]. On 
the other hand, the fact that another combination of the 
branching ratios, 



n _ f Br[B° 
R = 



K+TT-] + Br[B° -> K-n+} \ t 



Br[B+ -> K°7T+] 
(0.89 ± 0.07) exp 



Br[B~ -> K°t:-} J r B o 



(3) 



the B factories; S<j,K s <Sj/^Ks- As we will show in this 
letter, the modified QCD penguin contributions would 
also help to deviate R c and R n in an indirect manner. 

In our computation, we apply the QCD factorization 
(QCDF) [15,16] which offers us an ability of estimating 
the hadronic matrix element of C?7 7 . We extend the pa- 
rameterization in [15] by including the SUSY contribu- 
tions. Then, Eq. (18) in [15] can be rewritten as 
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= P \e l 

+7T K L 



P[e i0p +e a e 
- e 3/2 ^(e- 17 - 

+ e 3/2 e^(e~^ - qe 



■vy 



A^ a 



qe q e 

-17 



q c e m «ce« 



(4) 
(5) 
(6) 
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(e" 17 - q c e i9 "a e iua )] (7) 



which satisfies the isospin relation: 



V2A^o^ - 



-A, 



■>TT°K- 



t+K-' 



The parameters 4> a ,(t>,4>T,u,uJc and 9p,6 q ,6 qc are the 
CP conserving (strong) and the CP violating phase, re- 
spectively. The parameters are written as: 



with -^B- = 1.086 ± 0.017 [11, is almost consistent to the 
SM prediction leads us to the so-called largely enhanced 
electroweak penguin (EWP) mechanism [8] among vari- 
ous New Physics (NP) scenarios. Possible NP contribu- 
tions to EWP have been studied (see e.g. in [9-11]) and 
also an impact on the future experiments of the K decays 
are investigated [8,12]. 

In this letter, we analyze the supersymmetric contri- 
butions to the B — > Ktt process in a model independent 
way using the mass insertion approximation. We will 
show that the Z penguin diagrams with chargino in the 
loop contribute to the EWP significantly for a light right 
handed stop mass. Furthermore, the subdominant color 
suppressed EWP can be also enhanced by the the electro- 
magnetic penguin (O77) with chargino in the loop. The 
gluino contributions modify mainly the chromomagnetic 
penguins (0 Sg ), i.e. the QCD penguins. As shown in 
[13,14], this modified QCD penguin contribution would 
be the key to understand another hint of NP discovered in 
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where A„/A c = e KM e n and Rr-k = ~ 1.01. 

Comparing to [15], we have three extra CP violating 
phases, 9p,9 q ,9 qc which may be induced by SUSY. 

contains the QCD (EW) penguin conto- 



rt 



4(3. EW;4,EW) 

bution of both SM and SUSY. The SUSY contributions 



can be obtained by replacing SM Wilson coefficients in 
to the SUSY Wilson coefficients. a\ and 



a 4,(3, EW;4,EW) 

a-i are the color allowed and color suppressed tree con- 
tributions, which do not contain SUSY contributions. 0? 
represents the so-called weak annihilation contributions. 
We use the default values of all the input parameters in 
[15] with pa = Ph = in our numerical analysis. As a 
result, the SM contributions within QCDF are given by 



(Pe i5p )*M 
(e 3/2 e^)sM 

(ge^) SM 



-0.0989e-* ao22 ,(e a e^) SM 
0.231 e - ia077 , (e T e^) S M 
0.136e-' 075 

, (gee ) sm - 

(e 3 / 2 e i0 )sM 



= -0.0202e l025 
= 0.216e* 00077 
0.0143e- i088 
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Notice that all the CP violating phases apart from 7 are 
zero in SM. These theoretical values lead to 

R c = 1.09(1.60), R n = 1.11(1.73), R = 0.96(1.46) (8) 

for 7 = 7r/3(27r/3). 

We first introduce a parameterization in which the 
SUSY contribution manifests itself. By assuming the 
same strong phases for SM and SUSY, which is not a 
bad assumption in QCDF where the strong phases enter 
as higher-order correction, we can write 
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' a 1 EW^ SVSY 
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\RKir(%3 EW ^4 Bff)sM 



l 4.EW 
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(9) 
(10) 

(11) 



(12) 
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The index SM (SUSY) means to keep only SM (SUSY) 
Wilson coefficients in EW ^ and /??/ EW y £a, £3/2 and 
€t also include the QCD and EW penguin with the u in- 
dex. However, these are always suppressed by the factor 
€km — 0.020 comparing to the SUSY contributions with 
the c index. Therefore, neglecting this small contribu- 
tion, the SUSY effects modify £{ a ,3/2, t} as 



(e a e ) 1 ^ ; , , e , , , (£3/26 v ) 



(e 3 / 2 e^) £ 



(e T e l4,T ) = 



|l + fce l ^| 
(ere^) s , 
\l + ke ie 'p 



(15) 



In order to have a general picture of the B — > Ktt 
puzzle, an expanded formulae in terms of e-r, £bw = 
q x e 3 / 2 and e EW = qc x er is often useful. By assuming; 
i) the strong phases are negligible, i.e., (j) a ,(f>,ui,ujc are 
all zero, ii)the annihilation tree contribution is negligible, 
i.e. e a ~ 0, iii)the color suppressed tree contribution is 
negligible, i.e. e 3 /2e*^= eye^we can write R c and R n as 
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FIG. 1. Scattered plot for R n versus R n with 
{k = 0,1 = l,ra = 0} (circle (blue)), {k = 0, 1 = 2, m = 0} 
(triangle (red)) and {k = 0.3, / = 1.5, m = 0} (square (green)) 
by varying j,0',9'„. The ellipse indicates the experimental 
value. 
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(16) 
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Now, let us find the configuration which leads to 
Rc — Rn ft 0.2. From Eq. (17), we can find that in gen- 
eral, the larger the values of ex, £ew and £% w are, the 
larger the splitting between R c and R n we would acquire. 
Considering the SM value of qc xej, the second term in 
Eq. (17) has only a tiny impact unless m is of order 
10. We can see that the phase combinations Op — q and 
6p + 7 also play an important role. The possible solu- 
tion to the R c — R n puzzle by enhancing 6ew, which we 
have parameterized as I, has been intensively studied in 
the literature [6-8]. As we will see in the following, er 
can also be enhanced when ke lBp term contributes de- 
structively against the SM and diminish P (see Eq. (9)). 
However, since P is the dominant contribution to the 
B — > Ktt process, the branching ratio is very sensitive to 
ke lSp . Therefore, we are allowed to vary ke l0 p only in 
the range of the theoretical uncertainty of QCDF, which 
gives about the right sizes of the B — > Ktt branching ra- 
tios. Hence, we would be able to reduce P at most by 30 
%, which can be easily compensated by the error e.g. in 
the transition form factor F B ^ n ' K . 

Considering the tiny effect from the second term in 
Eq. (17), in order to achieve R c — R n > 0.2, we need 
£t£ew larger than about 0.1 or equivalently, eew larger 
than about 0.5 with e^ M . In figure 1, we show scattered 
plots by varying j,0' p ,0' q in a range of n/3 < 7 < 27r/3 
and —it < 0' < 7r (interval of 0.2 for each phase). 

4(C) 

It should be noted that the figure is obtained by us- 
ing the full formulae of R c and R n with the full QCDF 
result. The circles (blue), triangles (red) and squares 
(green) represent the result with {k = 0, 1 = l,m = 0}, 
{k = 0,1 = 2, to = 0} and {k = 0.3,1 = 1.5, to = 0}, 
respectively. We can see that with k = 0, we need / > 2 
to reproduce the experimental values while an inclusion 
of a small amount of k lowers this bound significantly. 



Now we investigate the SUSY contribution to the 
B — > Kir and the possible values of the parameters k 
and I to be reached. We work in the mass insertion ap- 
proximation which is a technique to evaluate the SUSY 
contributions in a model independent manner. In this ap- 
proximation, the extra FCNC as well as the new sources 
of CP violation enter in terms of the complex parameters 
(dfj)AB called mass insertion (MI) where A, B are left (L) 
or right (R) indices, i, j are generation indices and q refers 
the up or down sector. While we will here try to con- 
strain the MI from the experiments, it is also computed 
by given soft SUSY breaking models. As mentioned in 
the introduction, there are two kinds of new diagrams, 
one from the gluino loop and the other from the chargino 
loop. As can be found in [17], the gluino contributions 
to the Wilson coefficients are proportional to {& 23 )ll(rr) 
apart from the magnetic terms, Cj-y and C$ g , which ad- 
ditionally receive contributions proportional to the chi- 
rality flipping {8 23 )lr(rl)- For the moderate SUSY mass 
configurations we will consider in the following, we find 
that only the (^23) lr(rl) terms, which is enhanced by a 
factor nig/mb comparing to the (^23) ll(rR) terms, play 
a significant role. As the chargino contributions, we con- 
sider the general SUSY models with light chargino mass 
as well as the case with light right handed stop. As given 
in [14], the Wilson coefficients of the chargino diagrams 
are proportional to the up sector MI and the leading term 
come from {S 32 ) RL and {S% 2 ) LL . The others {S^l^rr), 
(^32) RL and (£32) rr are Cabbibo suppressed by order 
A, A 2 , A 3 , respectively where A ~ 0.22. The general for- 
mulae for the chargino Wilson coefficients are given as 



- (^32) liAf L 



+ (<&) 



RL 



\r T)RL 



(18) 



where the loop functions R F L and Rp L can be found in 
Rcf. [14]. The index F refers to M 1 (electromagnetic- 
penguin ), Af 9 (chromomagnetic-penguin), and C (Z- 
pcnguin ). It is important to note that the functions 
Rfrrr.g depend on the bottom Yukawa coupling so that 
it can be enhanced for large tan (3. Also the function R^j L 
largely increases when decreasing the mass of the right- 
stop. In the following, we consider only dominant three 
mass insertions, (5 23 )lr, {$32)ll and (5 32 )rl discussed 
above, which have some enhancement factor and have a 
potential to lead to a large SUSY contributions to the 
B — > Kit process. 

Now let us show our numerical result. As a default 
value of the SUSY parameters, we chose 

m g = 500GcV, rhg = 500GeV, 
m fjJ = 125GcV, m 2 = 150GeV, ix = 250GeV (19) 

First we present SUSY contributions to each a? ^ EW ^ 

from each MI. The (S 23 )lr and (5 32 )ll terms which con- 
tribute to Cr 7 and C% g lead to 



c, o 



-36.7(<5&) Lfl , 



c ,9 
y A,EW 

C,SM 

a 4,EW 



27.7(6* 3 ) LR , 



a 



C,SM - 



-0.00110 tan f3(5% 2 ) LL 



a 



4,EW r 

1 C,SM " 

a A.EW 



;0.148tan/?(«J£ 2 ) LL , 



respectively and the (5 32 )rl term which contributes to 



Z penguin are obtained as 

a c ' z a c ' z 

-4^ ~ 1.68(S$ 2 )hl, 4^ * 1.18(%)«l. 

a 4,EW a 3,EW 

Collecting all the SUSY contributions, our parameters 
k, I, and m are obtained as 

ke l9p = -0.0019 tan/3(%) iL - 35.0(<5^) ifl + 0.061 (^^ 
le"* = 0.0528 tan (3{5^ 2 ) LL - 2.78(6$ 3 ) L r + 1.11(<5&)jil 
me ie «c = 0.134 tan/3(5&) iL + 26.4(^ 3 ) Lfl + 1.62{5% 2 ) RL 

Note that we do not consider (S 23 )rl here but it is the 
same as (S 23 )lr with an opposite sign (see also [18]). 

Let us first discuss the contributions from a single 
mass insertion (S^ 2 ) LL , {5 23 ) L r or {5% 2 ) RL to {k,l,m}; 
keeping only one mass insertion and switching off the 
other two. Note that, as is well known, the absolute 
values of mass insertions (S 32 )ll and (S 23 )lr receive 
constraints from the b — > 57 branching ratio which are 
I tan/3 x {6% 2 ) LL \ < 1 and \(5$ 3 ) LR \ < 0.005. The mass 
insertion {S 32 )rl has only a constraint from its defini- 
tion K^)^! < 1- Firstly we discuss the (S 32 )rl term. 
Using K^)/?!,! = 1, the maximum value is found to be 
{k,l,m} = {0.061,1.11,1.62}. Thus, in this case where 
k is almost negligible, we would need I ~ 2 to explain 
the experimental data. We have a chance to enlarge the 
coefficients for (5 32 )rl by, for instance, increasing the 
averaged squark mass rhq. However, even if we choose 
fhq = 5 TeV, we find that I is increased only by 20 to 
30 %. Secondly, let us evaluate the (S 23 )lr and (6 32 )ll 
terms. Including the constraints from the b — > S7 branch- 
ing ratio, the maximum contributions from {5 23 )lr and 
(5 32 )ll are found to be {k,l,m} = {0.18,0.014,0.13} 
and {0.0019,0.053,0.13}, which are far too small to ex- 
plain the experimental data. The coefficients for (S 23 )lr 
depend on the overall factor 1 /fhq and on also the vari- 
able of the loop function x = m g /fhq and we found that 
nig = rhq = 250 GeV can lead to 100 % increase. How- 
ever, the value of I is still too small to deviate R c — R n 
significantly. As a whole, we found that it is extremely 
difficult to have R c — R n > 0.2 from a single mass inser- 
tion contribution. 

Let us try to combine two main contributions, (S 23 )lr 
and (S 32 )rl terms. Using the default SUSY masses 
in Eq. (19) and including the b — > sj constraint to 
I (^23)^^1 ) the maximum value is found to be {k, l,m} — 
{0.24, 1.12, 1.48}. The resulting R c — R n are given as cir- 
cles (blue) of Fig. 2 by varying arg(<5 2 i 3 )Lfl, arg(<^ 2 ) flL , 7. 
We can see that the experimental data are not repro- 
duced very well. As discussed above, for a large value 
of the averaged squark masses, I increases while k de- 
creases. On the contrary, k also depends on the ratio 
of gluino and squark masses. Hence we need to optimize 
these masses so as to increase k and I simultaneously. For 
instance, with nig =250 GeV and fhq=\ TeV, we obtain 
{k, I, m}={0.30, 1.36, 1.90} with which we find that quite 
a few points become well within the experimental bounds 
of R c and R n as shown in triangles (red) of Fig. 2. 

Before concluding, let us mention about the relaxation 
of the b — > S7 constraint as a possible solution when the 
experimental values of R c and R n remain around their 



FIG. 




= 1 



2. Results of SUSY models with \(5% 2 ) RL \ 
and |(523)-L-r| = 0.005 (circle (blue) and triangles (red); 
see the text for the other parameters) and |(5aa)-Ri| = 1, 
\($23)lr\ = 0.018 and |tan/3 x (5% 2 ) LL \ = 6 (square (green)). 

current central values. The branching ratio of b — > 
depends not only on the absolute value of the mass in- 
sertion but also on the argument of it due to the over- 
lap term of the SM and SUSY contributions. The con- 
straints to the absolute value considered above are ob- 
tained independently from the phase. On the contrary, 
for a certain value of the argument, much larger absolute 
value could be allowed. For example, | (^23) = 0-02 
with arg(<5 23 )ift — ±tt and |tan/3 x (5 32 )ll| = 6 with 
arg((5g 2 )LL — can also accommodate the experimen- 
tal value of the 6—5-37 branching ratio. However, 
this does not help unfortunately since the constraint to 
l(^32)ii| is s tiH t 00 strict to deviate R c and R n signifi- 
cantly and also the phase condition for the (^23)^^ term, 
arg((5 23 )£fl — i 71 ") works in a direction of diminishing 
Rc — Rn- The situation is different once we consider 
the (623) lr an d (£32 )ll terms simultaneously. When a 
condition | (^23) z-Jel = O(10 2 ) x | tan/3 x (8% 2 )ll\ is satis- 
fied, the cancellation between the (£| 3 ) lr an d (^ 3 2)ll 
terms occurs. Under this circumstance, we find that 
the same kind of relaxation for the absolute values can 
be expected but for much larger range of the phases 
comparing to the previous case. The squares (green) 
in Fig. 2 show a scattered plot with HS^rl] = 1, 

IO&UhI = °- 018 and I tan/ 3 x ( S 32)ll\ = 6 by vary- 
ing the phases in a rage of 7r/3 < 7 < 2ir/3 and 
-7T < arg(5 3 l 2 ) fl ,L,arg(5 2 z 3 ) Li? ,arg((5 3 l 2 ) LL < n and ex- 
cluding the points which do not satisfy the b —> sj branch- 
ing ratio. We observe much larger J? c — R n in this scenario. 

In conclusion, we have examined the supersymmet- 
ric models as a solutions to the B — » Kir puzzle. We 
have shown that the Z penguin diagram with chargino 
in the loop could be enhanced for the small value of 
right handed stop mass and order one (£32) rl- We, 
however, found that this contribution itself is not large 
enough to solve the puzzle when choosing moderate val- 
ues of the SUSY particle masses. We also found that 
B — > Kit receives a large gluino chromomagnetic pen- 
guin contributions ((S^lr term), which could explain 
another hint of new physics discovered in the B facto- 
ries, S^Ks < Sj/^ Ks . While the (6$ 3 ) lr term mod- 



ifies only the QCD penguins and does not solve the 
Rc — Rn puzzle directly, we found that it plays a com- 
plementary role to the (S^RL term and we can well 
explain the experimental values of R c and R n within 
the b — ► S7 constraint. The chargino electromagnetic 
penguin diagram could also enhance EWP especially for 
large value of tan/3, however, we found that the stringent 
constraint to | tan/3 x (<5 32 )ll| from the branching ratio 
of b — > sj prevents the (<5 32 ) ll term influencing the val- 
ues of R c and R n . We examined a possible relaxation of 
the b — > sj constraint considering the phase configura- 
tion of the mass insertions. We found that the relaxation 
is quite possible especially, when a cancellation between 
the {823) lr and (<5 32 )ll term occur. 

As represented in the relaxation scenario, we found 
that the R c and R n data and furthermore the full spec- 
trum of the B — > Ktt branching ratios constrain the 
SUSY parameters for FCNC and CP violation very 
severely, which would provide a great opportunity of im- 
proving our knowledge of the SUSY breaking nature [19]. 
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